First results on angular distributions of thermal dileptons in nuclear collisions 
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The NA60 experiment at the CERN SPS has studied dimuon production in 158A GeV In-In 
collisions. The strong excess of pairs above the known sources found in the complete mass region 
0.2<M<2.6 GeV has previously been interpreted as thermal radiation. We now present first results 
on the associated angular distributions. Using the Collins-Soper reference frame, the structure 
function parameters A, fi and v are measured to be zero, and the projected distributions in polar 
and azimuth angles are found to be uniform. The absence of any polarization is consistent with the 
interpretation of the excess dimuons as thermal radiation from a randomized system. 



Lepton pairs are a particularly attractive observable 
to study the hot and dense matter created in high- 
energy nuclear collisions. Their continuous emission, 
undisturbed by final-state interactions, probes the entire 
space-time evolution of the fireball, including the early 
phases with the conjectured QCD phase transitions of 
chiral symmetry restoration and parton deconfinement. 
To the extent that the bulk constituents of the expand- 
ing matter (hadrons and partons) equilibrate, the direct 
lepton pairs generated by them are commonly referred to 
as 'thermal radiation'. Our previous work has indicated 
'thermal' dilepton production to be largely mediated for 
M<1 GeV by ir + ir~ annihilation via the strongly broad- 
ened vector meson p jlj], and for M>1 GeV by partonic 
processes like qq annihilation 0, E|. The two dilepton 
variables basically explored in this work were mass M 
and transverse momentum p-r, where the correlations be- 
tween the two were decisive in bearing out the nature of 
the emission sources in the two mass regions 0, 13] ■ 

Further information on the production mechanism and 
the distribution of the annihilating particles, comple- 
mentary to that from M and p<r, can be obtained from 
the study of dilepton angular distributions. This Letter 
presents the first measurement of full dilepton angular 
distributions in the field of high-energy nuclear collisions. 
It is restricted to the mass region M<1 GeV, due to the 
lack of sufficient statistics for M>1 GeV. The question 
asked is simple: can one get direct experimental insight 



into whether the radiating matter is thermalized? 

Historically, the interest in angular distributions of 
continuum lepton pairs was mostly motivated by the 
study of the Drell-Yan mechanism, following in particular 
the insight that the 'naive' QED interpretation Jl had to 
be significantly modified due to QCD effects [I Jo, 0, 0] • 
The differential decay angular distribution in the rest 
frame of the virtual photon with respect to a suitably 
chosen set of axes, ignoring the rest mass of the leptons, 
can quite generally be written as 

oc (1 + \cos 2 9 + /j,sin29coscj) H — sin 2 9cos2(/)) (1) 

a ail 2 

The angular dependence on polar angle 9 and azimuth 
angle 4> dates back to but the specific coefficients A, 
fi and v, the 'structure function' parameters, follow the 
nomenclature in e.g. [EES]- They are directly related to 
the helicity structure functions Wi defined in and in 
particular to the spin density matrix R7. of the virtual 



photon 7*, the main object of the spin analysis 0, HI 111- 
We have chosen here the Collins-Soper (CS) reference 
frame [f| , where the quantization axis z is defined as the 
bisector between the beam and negative target momenta 
Pbeam and —ptarget, which define the reaction plane. The 
polar angle 9 is then the angle between the momentum 
of the positive muon p^+ and the z axis, which define 
the decay plane, while the azimuth angle (f> is the angle 
between the reaction and the decay planes. However, the 
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FIG. 1: Isolation of an excess above the electromagnetic de- 
cays of neutral mesons (see text). Total data (closed circles), 
individual cocktail sources (solid lines) , difference data (thick 
triangles), sum of cocktail sources and difference data (line 
through the closed circles). Open charm still contained. 



are stable with respect to the bin widths chosen. 

The assessment of the two background sources and 
open charm is extensively discussed in the ref. [IJ. 
The centrality-integrated net mass spectrum after back- 
ground subtraction is shown in Fig. Q] together with the 
contributions from neutral meson decays: the 2-body de- 
cays of the rj, u> and c/> resonances, and the Dalitz decays 
of the r), rj and u. The data clearly exceed the sum 
of the decay sources. The excess dimuons are isolated 
by subtracting them from the total (except for the p), 
based solely on local criteria [HQ. The excess for M<1 
GeV is interpreted as the strongly broadened p which 
is continuously regenerated by 7r + 7r~ annihilation [H, @]. 
Two adjacent mass windows indicated in Fig. [T] are used 
for the subsequent angular distribution analysis: the p- 
like region 0.6<M<0.9, and the low-mass tail 0.4<M<0.6 
GeV. To exclude the region of the low-mx rise seen for all 
masses \2\, a transverse momentum cut of pt>0.6 GeV 
is applied, leaving about 54000 excess pairs in the two 
mass windows. The subtracted data for the uj and </>, 
about 73 000 pairs, are subject to the same further anal- 
ysis steps as the excess data and are used for comparison. 



particular choice of the frame is not really relevant here, 
since the determination of the full set of coefficients A, p 
and v allows to compute them in any other frame by a 
simple transformation This would not apply if only 
the cosQ distribution would be measured. 

In principle, the dilepton angular distributions can 
be anisotropic, with all structure function parameters 
A, /i, v ^ 0, ID, 0, H, EH- Even for spinless particles 
in the initial state like in 7r + 7r _ annihilation, the param- 
eters can still have any value ^ 0, since the spin density 
matrix of the virtual photon also receives contributions 
from orbital angular momentum [12 ]. Very elementary 
examples are qq and ir + ir~ annihilation along the beam 
direction for p-r=0. Here p—v—0 and A=+l for qq (like 
lowest order DY and A=— 1 for 7r + 7r~ [jjj], corre- 
sponding to transverse and longitudinal polarization of 
the virtual photon, respectively. However, a completely 
random orientation of annihilating partons or pions in 3 
dimensions (but not in 2 [ill ]) would lead to \,p,v = 
11 . T2, EH i and that is the case of prime interest here. 







Details of the NA60 apparatus are contained in 0, [3] • 
The data sample for 158 A GeV In-In collisions is the 
same as used in [H, 0, and the different analysis steps 
follow the same sequence: assessment of the combinato- 
rial background from tt and K decays by a mixed-event 
technique, assessment of the fake-matches (associations 
of muons to non-muon vertex tracks in the Si pixel tele- 
scope), isolation of the dimuon excess by subtraction of 
the known meson decay sources and charm from the net 
opposite-sign sample, and finally correction for accep- 
tance and pair efficiency. All steps are now done inde- 
pendently in each [ g^ggz ]-tj bin. The binning is varied 
depending on the goal, thereby assuring that the results 



The correction for the acceptance of the NA60 appa- 
ratus requires, in principle, a 5-dimensional grid in M- 
PT-y-cos8-cj> space. To avoid large statistical errors in 
low-acceptance bins, the correction is performed in 2- 
dim. cos9-(f> space, using the measured data for M, pr Q 
and y [l5| as an input to the Monte Carlo (MC) simu- 
lation of the cos9-(f) acceptance matrix. The sensitivity 
of the final results to variations of this M-p<r-y input has 
been checked, and the effects are found to be consider- 
ably smaller than the statistical errors of the results. The 
MC simulations were done in an overlay mode with real 
data to include the effects of pair reconstruction efficien- 
cies. The product acceptance x efficiency is illustrated in 
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FIG. 2: Spectrometer acceptance as a function of the two 
angular variables \cos6\ and |c/>|. 
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FIG. 3: Polar angle distributions of excess dileptons and of the vector mesons u> and 



Fig.dfor 0.6<M<0.9 GeV and p T >0.6 GeV. The rapid- 
ity coverage is 3.2<y<4.2 (+0.3<y cm <+1.3). 

The results on the angular distributions have been 
analysed in three different ways, distinguished by the 
method and the associated statistical/systematic errors. 
In the first and most rigorous method (1), the 3 struc- 
ture function parameters A, fi and v are extracted 
from a simultaneous fit of the 2-dimensional data on 
the basis of Eq.{T]), using a 6x6 matrix in the range 
-0.6<cos6»<+0.6 (bin width 0.2) and -0.75<co S( /)<0.75 
(bin width 0.25). The restrictions in range are enforced 
by regions of very low acceptance in the 2-dimensional 
acceptance matrix, masked in the projections of Fig. [2] 
The fit values are summarized in Table U for all 4 cases, 
the two excess mass windows, the u> and the cf>. Within 
errors, they are all compatible with zero. It is reassur- 
ing to see that this is also true for [i, as expected for a 
symmetric collision system at midrapidity on the basis of 
symmetry considerations [111 ]. 

In the second method (2), setting now /i=0, the 
2-dimensional acceptance-corrected data are projected 
onto either the |cos#| or the \§\ axis, summing over the 
two signs. The polar angular distribution is obtained by 
integrating Eq.([T]) over the azimuth angle (<j>) 

^L=A (l + \co S 2 9), (2) 
acosa 

while the azimuth angular distribution is obtained by 
integration over the polar angle (cos9) 

dN 



Ai{l + ^ + ^cos20) 



(3) 



The structure function parameters A and v can then be 
determined independently by I-dimensional fits to the re- 
spective projections. The data of the polar angular dis- 
tributions together with the fit lines according to Eq. {2J 
are shown in Fig. [3] for all four cases, using now 8 bins in 
\cos9\ (bin width 0.1). The distributions are seen to be 
uniform, and the fit parameters A, included in Table [H 
are again compatible with zero, within errors. To deter- 
mine the parameter v, A (contained in Eq.f3|)) is set to 
the measured value of A ± ax ■ The fit results for v on the 
basis of Eq. (J3|) , keeping the small number of bins used in 
method 1, are again zero, within errors (see Table [TJ> . 

In the third method (3), the inclusive measured dis- 
tributions in cosQ and <f> are analysed. A 1-dimensional 
acceptance correction is applied in each case, determined 
by using (as now measured) uniform distributions in <p 
(for cosQ) and in cosQ (for 4>) as an input to the MC sim- 
ulations. The number of bins in |cos#| is kept, while that 
in \<j)\ is increased to 10 (bin width 0.3). The data for 
the azimuth angular distributions together with the fit 
lines according to Eq.Q are shown in Fig. 2] The dis- 
tributions are again uniform, as are those for cos9 (not 
shown, since hardly distinguishable from Fig. [3]) . The re- 
sulting fit parameters for A and v (taking account again 
of A ± ax) are included in Table Q] As expected, the er- 
rors are smaller than for the other two methods, but the 
values of A and v are still close to zero, within errors. 

Figs. [3] and H] also contain the systematic errors at- 
tached to the individual data points. They mainly arise 
from two sources. The subtraction of the combinatorial 
background, with relative uncertainties of 1% [H, 0, 0], 
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FIG. 4: Azimuth angle distributions of excess dileptons and of the vector mesons ui and 



leads to errors of 2-3% of the net data for the kinematic 
selection used here. The subtraction of the meson de- 
cay sources causes (correlated) errors for the excess and 
the vector mesons u> and cf>. With respect to the excess, 
they range from 4-6% up to 10-15% in the low-populated 



TABLE I: Summary of results for pt>0.6 GeV on the struc- 
ture function parameters A, fi and v in the CS frame, ex- 
tracted from three different methods (see text). The \ 2 /ndf 
of the fits varies between 0.8 and 1.2. For a cut pr>1.0 GeV, 
the results are the same, within errors. 
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cos6-<f> matrix bins. This variation is well visible for the 
overall errors plotted in Figs. [3] and [4] Assuming, very 
conservatively, these errors to be uncorrelated from point 
to point, the (statistical) fit errors for A and v quoted in 
Table U would increase by 15-20% if the systematic errors 
would be added in quadrature. Further confidence into 
the stability of the results is obtained from their inde- 
pendence of the methods and the bin widths used. 

The global outcome from our analysis of dilepton an- 
gular distributions is straightforward: the structure func- 
tion parameters A, [i and v are all zero within the statisti- 
cal and systematic errors, and the projected distributions 
in polar and azimuth angle are all uniform. This applies 
not only for the excess dileptons as anticipated if of ther- 
mal origin, but also for the vector mesons ui and <p. While 
there may be a rather direct connection between the two 
findings in nuclear collisions, it is of interest to note that 
the result of A=0 has been reported before for p and ui 
production in p-p [l6| and ir~-C [l7j ]. 

We conclude, following the primary motivation of this 
study, that the absence of any polarization is fully con- 
sistent with the interpretation of the observed excess 
dimuons as thermal radiation from a randomized system. 
While this is a necessary condition, it is not sufficient. 
However, together with other features like the Planck-like 
shape of the excess mass spectra [!, UH , the exponential 
shape of the spectra EH and the global agreement 
with theoretical models both as to spectral shapes and 
absolute yields [D, Ell , the thermal interpretation has be- 
come more plausible than ever before. 
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